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We report on a study of cold reactive collisions between sympathetically-cooled molecular ions and 
laser-cooled atoms in an ion-atom hybrid trap. Chemical reactions were studied at average collision en- 
ergies (£co1i)Ab > 20 mK, about two orders of magnitude lower than has been achieved in previous 
experiments with molecular ions. Choosing N^+Rb as a prototypical system, we find that the reaction 
rate is independent of the collision energy, but strongly dependent on the internal state of Rb. Highly effi- 
cient charge exchange about four times faster than the Langevin rate was observed with Rb in the excited 
(5p) 2 P3/2 state. This observation is rationalized in terms of a capture process dominated by the charge- 
quadrupole interaction and a near resonance between the entrance and exit channels of the reaction. Our 
results provide a test of classical models for reactions of molecular ions at the lowest energies reached 
thus far. 



The combination of radiofrequency (RF) ion traps with 
magneto-optical or optical-dipole traps for the simultane- 
ous confinement of cold ions and atoms has recently en- 
abled the study of ion-neutral interactions at temperatures 
as low as a few millikelvin |Q]-0. Hybrid trapping tech- 
niques have opened up perspectives for the exploration of 
new mesoscopic quantum systems [H, for quantum inter- 
faces between atoms and ions [9] and for new methods 
to cool ions to ultralow temperatures [6][10]]. Moreover, 
chemical experiments with ions El 13 |7J [Til have started 
to approach an energy regime so far restricted to neutrals 
H2HT51 in which the quantum character of the collision can 
strongly influence the chemical dynamics llT6] - [T8ll . 

Until now, studies of reactive processes in hybrid traps 
have been restricted to atomic ions and neutral atoms. De- 
spite the apparent simplicity of these collision systems, a 
diverse range of chemical phenomena has been observed. 
The effects range from the prominent role of light-assisted 
processes in the radiative formation of molecules and the 
promotion of charge exchange in Ca + +Rb [4], the occur- 
rence of unusually fast reaction rates in Yb + +Ca and 
the influence of the ion's hyperfine state on the reactivity 
in Yb + +Rb Q. However, experiments with atomic species 
can only cover a fraction of the wealth of chemical phe- 
nomena, and it is necessary to extend studies to molecules 
to be able to fully explore the diversity of reactive effects 
at ultralow energies. 

In the present work, we extend ion-neutral hybrid- 
trapping techniques to enable the study of cold reactive 
collisions with molecular species by sympathetically cool- 
ing molecular ions via their interaction with laser-cooled 
atomic ions fT9l and immersing them in a cloud of ultra- 
cold atoms confined in a magneto-optical trap (MOT). To 
our knowledge, this is the first time that chemical reac- 
tions with molecular ions have been studied at collision 
energies down to (E co x)/k^ « 20 mK, about two or- 
ders of magnitude lower than has been achieved in pre- 
vious cold-collision experiments ll20l - [22l . By choosing 
sympathetically-cooled ions and laser-cooled Rb atoms 
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FIG. 1 . Experimental false-color laser-cooling-fluorescence im- 
ages of Coulomb crystals recorded over the course of reaction 
with ultracold Rb atoms. Sympathetically-cooled ions are made 
visible in molecular dynamics simulations. Color code: blue: 
laser-cooled Ca + ions; yellow: reactant ions; green: Rb + 
product ions; magenta: CaH + , CaOH + . See text for details. 



as a prototypical system, we characterize important fea- 
tures of cold reactive collisions between molecular ions and 
alkali atoms. We observe charge exchange between the two 
species and a strong dependence of the reaction rate on the 
electronic state of Rb. Whereas reactions with Rb atoms 
in the (5p) 2 Pz/2 state occur fast consistent with a clas- 
sical capture model, reactions with ground-state Rb atoms 
were found to be at least two orders of magnitude slower. 
The fast excited-state rate is explained by a long-range cap- 
ture of the atom by the molecular ion dominated by the 
long-range charge-quadrupole interaction and a near reso- 
nance between the electronic entrance and exit channels of 
the system at short range. Our results highlight the impor- 
tance of resonant effects in cold ion-neutral collisions and 
provide a test of classical capture models for reactions of 
molecular ions at the lowest energies reached thus far. 

Our hybrid-trapping apparatus has already been de- 
scribed in Ref. J4). In the present study, ions were 
generated in a linear radiofrequency (RF) ion trap from 
room- temperature background nitrogen gas using [2+1] 
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resonance-enhanced multiphoton-ionization (REMPI) via 
the Q branch of the vibrationless X 1 S+ —> a" 1 S+ 
transition |[23ll . The ions were subsequently thermalized 
through collisions with the background gas ['23']. 

After their generation, the ions were sympatheti- 
cally cooled by the interaction with laser-cooled Ca + ions 
to form Coulomb crystals fT9l . Images of the crystals 
were obtained by collecting the spatially resolved laser- 
cooling fluorescence of the Ca + ions using a CCD cam- 
era coupled to a microscope. The number of laser- and 
sympathetically-cooled ions in the crystals were deter- 
mined by comparisons of the experimental images with 
molecular-dynamics (MD) simulations as described in 
Refs. 12411251 . 

After ion loading, the Coulomb crystals were translated 
along the axis of the linear RF trap using a static electric 
field to be overlapped with the ultracold atoms. Clouds 
of ultracold Rb atoms at densities n R b = 4 x 10 8 cm -3 
were generated by laser-cooling in the MOT overlapped 
with the ion trap J4). Rate constants for reactive colli- 
sions between ions and atoms were determined by mon- 
itoring the decrease of the number 7V N + of reac- 
tant ions as a function of the reaction time t and fit- 
ting the data to an integrated pseudo-first-order rate law 
ln[iV N +(£)/iV N +(£=0)] = — fcpfot. Second-order rate con- 
stants k were determined from the pseudo-first-order con- 
stants fcpfo using k = &pfo/n R b. Both the sympathetically- 
cooled and the laser-cooled Ca + ions react with Rb. 
However, the total rate for Ca + +Rb is two orders of magni- 
tude smaller (4) than for N^+Rb so that the reorganization 
of the Coulomb crystals caused by loss of Ca + could be 
neglected over the timescale of the rate measurements for 
N+. 

Fig. [T](a)-(d) shows a sequence of false-color Coulomb 
crystal images and their MD simulations recorded over the 
course of a typical reaction experiment. The laser-cooling 
fluorescence of the Rb atoms was blocked by a color fil- 
ter so that only the Ca + fluorescence (blue) is visible in 
the images. The ions can only indirectly be seen as a 
dark non- fluorescing region in the center of the crystals, but 
have been made visible in yellow in the simulations. Ini- 
tially, Ca + Coulomb crystals containing ^1000 ions were 
loaded with typically ^25 sympathetically-cooled ions 
(Fig. [l](a)). Over the course of the experiment, the Njj~ ions 
were removed from the crystals by charge-exchange colli- 
sions with Rb atoms resulting in a reduction of the dark 
crystal core (Figs. [T](b)-(d)). The Rb + product ions (green 
in the simulations) remained trapped and were sympatheti- 
cally cooled by the Ca + ions. The magenta-colored ions in 
the simulated images represent CaH + and CaOH + which 
are impurities formed by reactions of Ca + with background 
H 2 and H 2 gas, respectively. 

To prove the chemical identity of the reaction products, 
resonant-excitation mass spectra of the Coulomb crystals 
were recorded [26]. Briefly, an additional RF drive voltage 
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FIG. 2. Resonant-excitation mass spectra of Coulomb crystals 
(insets) loaded with sympathetically-cooled ions (i) before 
reaction with ultracold Rb atoms and (ii) after reaction of the 
ions with Rb. The dashed vertical lines indicate the single-ion 
resonance frequencies of the relevant ion species. The weak fea- 
ture around 70 kHz in (ii) is identified with the Rb + product ions 
generated by charge exchange of with Rb. (iii) Spectrum of a 
crystal prepared without ions recorded after the same time of 
reaction with Rb as in (ii) illustrating the effect of the slow back- 
ground reaction Ca + +Rb— ^Ca+Rb + . See text for explanation. 



applied to the trap electrodes was used to resonantly ex- 
cite the radial motion of specific ion species. The motion 
of the excited ions couples to the laser-cooled Ca + ions 
via the Coulomb interaction and leads to a modulation of 
their laser-cooling fluorescence around the relevant reso- 
nant excitation frequencies. Fig. [2] shows mass spectra of 
representative Coulomb crystals (i) before and (ii) after re- 
action. Additionally, a pure Ca + crystal (iii) was reacted 
for the same duration as in (ii) to illustrate the effect of 
the slow background reaction Ca + + Rb. The insets show 
images of the corresponding crystals and the vertical lines 
indicate the theoretical single-ion resonance frequencies of 
the relevant ion species. 

Whereas the spectrum recorded before reaction only 
shows a single resonance corresponding to the excitation 
of the Ca + ions, the spectrum after reaction in (ii) shows 
an additional weak feature around 70 kHz assigned to 
the sympathetically-cooled Rb + product ions. Note that 
the positions of the resonances in a Coulomb crystal are 
slightly shifted from the single-ion frequencies because 
of Coulomb-coupling effects [26]. Resonances for the 
CaOH + , CaH + and the reactant ions, as well as for the 
Rb + ions formed by the background reaction Ca + + Rb in 
(iii) could not be observed in the spectra under the present 
conditions. MD simulations revealed that their numbers 
were too small to exert any noticeable back action on the 
Ca + fluorescence at the excitation amplitudes used. Higher 
amplitudes would have led to the melting of the crystals 
around the position of the Ca + resonance and were there- 
fore precluded. 

Fig. [3] (a) shows the dependence of the rate constant 
on the collision energy. In the experiments, the Rb ki- 
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FIG. 3. (a) Rate constants for N^+Rb as a function of the col- 
lision energy. Each data point represents an average over three 
independent measurements, the error bars indicate the statistical 
uncertainty (2a). The insets show the images of the Coulomb 
crystals corresponding to the three data points indicated, (b) 
ion kinetic-energy distributions from MD simulations cor- 
responding to the data points with the lowest (black) and highest 
energies (grey) displayed in (a). The insets show the relevant 
simulated crystal images, (c) Rate constant as a function of the 
population in the Rb (5p) 2 ^3/2 state. 



netic energies ((E^ n ) /k# ~ 200/xK) were much smaller 
than the ion energies (> 20 mK) so that the collision en- 
ergies were entirely dominated by the contribution of the 
ions. The motion of the ions can be separated into two dif- 
ferent components, a slow secular (thermal) motion in the 
time-averaged pseudopotential of the trap and a fast micro- 
motion driven by the RF fields [19]. The secular energies 
((E sec ) /k B =12 mK for Ca+ and 14 mK for ) and micro- 
motion energies were determined from the ion trajectories 
obtained in the MD simulations of the crystal images l24l . 
The effective (total) kinetic energy E tot of the ions was gov- 
erned by their micromotion energy which vanishes for ions 
located on the central trap axis and increases quadratically 
with the distance lfT9l . The lowest achievable energies are 
limited by the precision with which the ions can be lo- 
calized on the axis. The axialization procedures adopted in 
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FIG. 4. Energy diagram of relevant electronic entrance and prod- 
uct channels for N^+Rb. The closest product channels associated 
with the N2"+Rb(5s) and (5p) entrance channels are indicated in 
green and red, respectively. 



the present experiments (to be described in detail in Ref. 
E3) allowed for a micromotion energy limit correspond- 
ing to 17 mK. 

For the measurements reported in Fig. [3] (a), Coulomb 
crystals with strings of nitrogen ions (like the ones shown 
in Fig. [T](a)) were aligned parallel to the trap axis so that 
all ions exhibit approximately the same micromotion en- 
ergy [24 j. The collision energies were varied by displac- 
ing the Coulomb crystals from the axis using static electric 
fields causing the ions to move into regions where they 
acquire an increased micromotion energy (see insets). The 
collision-energy distributions obtained from the MD simu- 
lations which correspond to the data points with the lowest 
and highest energies in Fig. [3] (a) are displayed in Fig. [3] 
(b) for reference along with the corresponding simulated 
images. As can be seen from Fig. [3] (a), within the uncer- 
tainty limits the rate constant is essentially independent of 
the collision energy in the interval studied. 

Because the ultracold atoms are constantly excited dur- 
ing laser cooling, a fraction of the reactions occurs with 
Rb atoms in the (5p) 2 Ps/2 state. The population in the 
excited state was adjusted by varying the intensity of the 
Rb cooling laser (4). The rate constant shows an ap- 
proximately linear increase with the (5^) 2 P3/2 popula- 
tion (Fig. [3](c)). In analogy to Ref. [4], a kinetic model 
of the form k — \ [k s (l + p s ) + k p p p ], where k s ^ p and 
p SiP denote the rate constants and populations in the Rb 
(5s, 5p) states, respectively, was fitted to the data. The fit 
yielded k p = 2.4(13) x 10 -8 cm 3 s _1 and an upper bound 
k s < 2 x 10 -10 cm 3 s _1 . The uncertainty of k p is domi- 
nated by a systematic error of « 50% in the Rb density. 

Because of the large difference between the ionization 
energies of N 2 and Rb, a minimum of 11.4 eV of energy 
are released in the reaction. By comparison, the rotational- 
vibrational energy of the ions (^ 0.025 eV for a room- 
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temperature distribution) as well as the collision energy 
(« 10 -6 eV) are small, so that the available energy for 
the reaction is entirely dominated by the electronic contri- 
bution. Fig. [4] shows an energy diagram depicting the rel- 
evant electronic entrance and product channels (see, e.g., 
Ref. [28 ] for a comprehensive overview of the electronic 
states of N 2 ). 

Previous studies using keV ion beams showed that 
charge exchange between alkali atoms and N J is most ef- 
ficient in a near-resonant transfer of the electron into a 
Rydberg state of the neutral product molecule which is 
built on an ion core with the same electronic configura- 
tion as the reactant ion ||29l . In the present case, the 
N2"(X + )+Rb(5s) entrance channel is closest in energy to 
a product channel forming N 2 in the C 3 T1 U electronic 
state (energy mismatch AE = 372 meV, see Fig. [4]). 
However, electron capture by N^ (X + ) (molecular orbital 
configuration (2a^) 2 (2a n ) 2 (l7r n ) 4 (3cr^) 1 ) to form N 2 in 
the C state ((2aJ 2 (2a n ) 1 (l7r n ) 4 (3aJ 2 (l7rJ 1 ) entails a 
significant re-arrangement of the electronic configuration 
of the molecular core. This process is therefore expected 
to be inefficient, in line with the present observation of a 
slow reaction rate in this channel. On the other hand, the 
excited (X + )+Rb(5p) entrance channel is near reso- 
nant with product channels forming N 2 in the the close- 
lying G 3 II n , D 3 £+ and F 3 II n states (energy mismatches 
A£"=183 meV, 151 meV and 8 meV). These heavily mixed 
states [ 30 ] are among the lowest Rydberg states of N 2 built 
on the X + 2 £+ and A + 2 Ii u ion cores so that capture 
of the electron into these states is expected to be highly ef- 
ficient, in agreement with the fast rate observed in the ex- 
cited channel. The Rydberg molecules subsequently pre- 
dissociate into the atomic fragments N( 4 5)+N( 2 Z)) lf29l 
which, however, cannot be detected in the present experi- 
ment. 

In previous studies on cold ion-neutral collisions in hy- 
brid traps 01 13, the rates were found to be in agreement 
with a limiting value set by Langevin theory 13T1 . The 
rate constant for the excited channel determined above 
(k p = 2.4(13) x 1CT 8 cm 3 s _1 ), however, is about four 
times larger than the Langevin value (k^ = 6.6 x 
10 -9 cm 3 s _1 ). This finding suggests that intermolecular 
forces other than the Langevin (charge-induced dipole) in- 
teraction play a role. In the present case, the most likely 
candidate is the interaction of the charge of the ion with the 
permanent electric quadrupole moment of Rb generated by 
the anisotropic charge distribution in the (5p) 2 Pz/2 state 

The rate constant was calculated with a classical 
model based on an interaction potential [33] incorpo- 
rating the Langevin and charge-quadrupole interactions. 
The quadrupole moment of Rb (5p) 2 P^/2 was computed 
within the single-electron approximation following Ref. 
Il34l . In line with the reaction mechanism discussed above, 
a capture approximation l35l was adopted, i.e., the reac- 
tion was assumed to be governed by the long-range inter- 



actions and to occur with unit probability once the reaction 
complex was formed. Within this framework, the classical 
reaction cross section = nb 2 was calculated from the 
maximum possible impact parameter b for which the height 
of the barrier in the centrifugally corrected potential does 
not exceed the collision energy 1551 . From a^ c \ the cor- 
responding rate constant at a collision energy E co \\/k^ — 
23 mK was obtained to be k^ = 1.7 x 1CT 8 cm 3 s"\ 
which agrees with the experimental value within its uncer- 
tainty limits. 

Notwithstanding the approximations inherent in the 
model (capture approximation, neglect of quantum effects 
and the molecular structure), the result underlines the ne- 
cessity to reach beyond the Langevin picture to explain the 
dynamics in the present system. It also suggests that - at 
least for the present case - classical capture models seem to 
adequately reproduce the kinetics in the millikelvin regime. 
This conclusion is consistent with the theoretical results 
of Refs. 010 [35j which predict a robust performance of 
classical capture models for ion-neutral reactions down to 
ultracold (<mK) temperatures for all but the lightest sys- 
tems. This behavior was rationalized in Ref. l35l in terms 
of a cancellation of tunneling and reflection effects at the 
centrifugal barrier. 

In conclusion, we have studied reactions between 
sympathetically-cooled ions and ultracold Rb atoms at 
average collision energies down to ^20 mK in an ion-atom 
hybrid trap. Our results highlight the importance of res- 
onant electronic effects in electron-transfer processes be- 
tween ultracold alkali atoms and cold molecular ions. A 
classical capture model was found to adequately describe 
the kinetics down to the lowest collision energies for chan- 
nels in which the short-range electron-transfer probability 
is expected to be close to unity. Whether short-range ef- 
fects such as the rotational- vibrational state of the molec- 
ular ion which are not accounted for in the capture picture 
play a role in the present system remains to be established 
in further studies performed with fully state- selected ions 
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